Electroweak phase transition in the MSSM with four generations 



S. W. Ham«, S. K. Oh^ 2 ), D. Son«, 

W Center for High Energy Physics, Kyungpook National University 

Daegu 702-701, Korea 
^ Department of Physics, Konkuk University, Seoul 143-701, Korea 



Abstract 

By assuming the existence of the sequential fourth generation to the minimal supersym- 
metric standard model (MSSM) , we study the possibility of a strongly first-order electroweak 
phase transition. We find that there is a parameter region of the MSSM where the elec- 
troweak phase transition is strongly first order. In that parameter region, the mass of the 
lighter scalar Higgs boson is calculated to be above the experimental lower bound, and the 
scalar quarks of the third and the fourth generations are heavier than the corresponding 
quarks. 
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I. Introduction 



As a mechanism to explain the baryon asymmetry of the universe, the electroweak baryogenesis is 
given wide attentions, since it can be tested in the future high energy experiments [1]. Sakharov 
have several decades ago established the three essential conditions for generating dynamically 
the baryon asymmetry of the universe from a baryon-symmetric universe [2]. As is well known, 
the three conditions are the presence of baryon number violation, the violation of both C and CP, 
and a departure from thermal equilibrium. The possibility of the electroweak phase transition 
has already been exhaustively studied, which can provide the baryon number violation and the 
violation of both C and CP. The remaining Sakharov condition, the departure from thermal 
equilibrium, may be fulfilled at a weak scale temperatures if the nature of the electroweak phase 
transition is first order. The difficulty of the standard model (SM) is that the strength of the first 
order electroweak phase transition, which must be strong enough for preserving the generated 
baryon asymmetry at the electroweak scale, appears too weak for the experimentally allowed 
mass of the SM scalar Higgs boson [3]. 

Thus, it seems that electroweak baryogenesis requires a new physics beyond the SM at weak 
scale [4]. The minimal super symmetric standard model (MSSM) has been studied intensively 
within the context of electroweak baryogenesis. It is observed that, if one of the scalar top quark 
has a mass smaller than top quark mass, the MSSM may possess a parameter region where the 
electroweak phase transition is strong enough [5]. In this scenario of a light scalar top quark, the 
requirement that the electroweak phase transition should be strongly first order is equivalent to 
the imposition of an upper bound of about 120 GeV on the lightest Higgs scalar boson mass of 
the MSSM. 

In any supersymmetric standard model, there are superpartners to ordinary quarks and 
leptons. Experimentally, no scalar quark or scalar lepton that is lighter than 5 GeV has been 
discovered. Thus, the ordinary quarks and leptons, up to bottom quark, are lighter than their 
corresponding superpartners. Possibly, top quark might be an exception, if a scalar top quark 
is lighter than top quark. 

The idea that a scalar top quark might be lighter than top quark is not in accord with 
the behavior of lighter fermions, but is allowed by present experiment and accommodated in 
the MSSM within the context of the electroweak phase transition. If the scalar top quarks are 
not degenerate in mass, the 2x2 mass matrix for the scalar top quarks yields a lighter scalar 
top quark and a heavier one in the MSSM. In order for the electroweak phase transition to be 
strongly first order, the mass of the scalar top quark should be either smaller than 170 GeV 
or larger than 1 TeV. Thus, the lighter scalar top quark should have a mass smaller than 170 
GeV. The heavier scalar top quark may have a mass comparable to the SUSY breaking scale, 
say, between 1 and 2 TeV. A light scalar top quark suggested in the MSSM scenario, which 
depends strongly on chargino and neutralino masses, might soon face experimental examination 
at Tevatron. 

The possibility of a scalar top quark heavier than top quark has been considered in scenarios 
other than the MSSM. For example, in the next-to-minimal supersymmetric standard model 
(NMSSM), the strongly first order electroweak phase transition takes place where the top quark 
mass is smaller than the mass of the scalar top quark. Within the framework of the MSSM, 
we need to enlarge the model somehow to accommodate a scalar top quark heavier than top 
quark. We examine the possibility for the strongly first-order electroweak phase transition in 
the MSSM by introducing an extra generation of fermions. 

In this paper, we study the effect of the fourth generation of quarks on the strength of the first 
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order electroweak phase transition in the MSSM. We find that introducing an extra generation 
of fermions might also lead to a scalar top quark heavier than top quark. Although the number 
of the SM neutrino species has already been fixed experimentally as three, one may find a lot of 
articles in the literature that mention about the fourth generation and study on the assumption 
of its existence [6,7]. In principle, the SM, as well as the MSSM, can accomodate any number 
of generations. In various contexts, the MSSM with four generations has been studied [8,9,10]. 

Our study shows that the fourth generation is found to enhance the strength of the first-order 
electroweak phase transition, while the scalar Higgs boson mass is calculated to be larger than 
the experimental lower bound, and the scalar quarks of the third generations are comparable to 
the supersymmetry breaking scale (Mstjsy = 1 TeV), in a reasonably wide region of parameter 
space in the MSSM. In our scenario, a light scalar quark is not necessarily required to ensure 
the first order electroweak phase transition be strong; the scalar quarks of the third and the 
fourth generations may be heavier than the corresponding quarks; and the scalar Higgs boson 
mass lies above the experimental lower bound. 

II. Higgs potential in decoupling limit without mixing 

Let us study a particular, yet reasonable as well as plausible, form of the Higgs potential in the 
MSSM with four generations of quarks for the electroweak phase transition. We consider only 
the third and the fourth generations, and assume that there is no mixing between them. The 
fourth generation appears simply in a repetitive manner. As is well known, there are two Higgs 
doublets in the Higgs sector of the MSSM, namely, Hf = (H$, H{) and H% = (H£,H§). After 
electroweak symmetry breaking emerge five physical Higgs bosons: two neutral scalar Higgs 
bosons (h ,H), one neutral pseudoscalar Higgs bosons (A), and a pair of charged Higgs bosons 
(H^. We assume that the mass of h is lighter than that of H. At the tree level, the Higgs sector 
of the MSSM depends on only two free parameters. We take them to be the ratio tan (3 = v-ilv\ 
of the two real vacuum expectation values (VEVs) v\ of H® and of and m^, the mass of 
A. In this paper, we assume that CP is conserved in the Higgs sector by choosing all parameters 
in the effective Higgs potential to be real. 

In the decoupling limit, where tua 3> rnz, with fixed tan/3, only one linear combination of 
the two neutral scalar Higgs bosons, 



remains light at the electroweak scale [11]. In this limit, the tree-level Higgs potential at zero 
temperature can be expressed in terms of <\> as 



Since all quartic terms have gauge coupling coefficients in the MSSM, the quartic Higgs self- 
coupling A is given as A = {g\ + g^)/^- Note that there is an upper bound on the mass of h 
as rrih < mz \ cos2/?| at the tree level in the MSSM. In this limit, the couplings of h to gauge 
bosons and fermions are identical to the couplings of the SM Higgs boson, which implies that 
one cannot distinguish phenomenologically the SM scalar Higgs boson from h [12]. Thus, one 
might expect that has the same experimental lower bound as the SM scalar Higgs boson in 
the decoupling limit [13]. The current experimental lower bound on the mass of the SM scalar 
Higgs boson is about 114.5 GeV. 



<f) = v^cos PRe(Hi) + v / 2sin/3Re(^) , 



(1) 



V o (cf>,0) = -m 2 cj> 2 + -cf> 4 . 



(2) 
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Now, at the one-loop level at zero temperature, the effective Higgs potential is given by 



V(0,O) = VoW>,O) + Vi(0,O) , 



(3) 



where the one- loop contribution V\{4>, 0) at zero temperature is obtained via the effective po- 
tential method as [14] 



nimf((j)) 
64tt 2 



log 



A 2 



(4) 



where I stands for various participating particles: the gauge bosons W, Z, the third generation 
quarks and scalar quarks t, 6, t±, I2, b\ and 62, as well as the fourth generation quarks and 
scalar quarks t', b f , t[, t' 2 , b[ and b' 2 . The renormalization scale in the above one-loop effective 
potential is set as A = mz- The degrees of freedom for each particle are: nw = 6, nz = 3, 



n t = n b 



-12, nr 



6 (i = 1,2 



dependent masses are given by m 2 ^ ((f)) 
m 2(0) = h 2 b cosp 2 (f) 2 /2, ml(4>) = 



n# = ny -- 
g 2 4> 2 /4, m 2 z (<P) 
h 2 ,smp 2 ^ 2 /2, m 2 b ,(<f>) 



12, Tip 



ny = 6 (i = 

~- (gl+g%)i 2 /4, m 2 (q 
h\, cos [3 2 (j) 2 /2, and 



1,2). Their field- 

i) = h 2 sin [3 2 <p 2 /2, 
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(5) 



with q = t,b, t', b' . In the above expression for the scalar quark masses [15], we have 



TO- 
IL 



= toq + to 2 (0) + ^- - - sin 2 cos 2/3to|( 
2 

= to 2 / + to 2 (</>) + - sin 2 9w cos 2/3to| (0) , 
| L (</>) = rn 2 Q + m\ (0) + (- X - + i sin 2 0^ cos 2/3m| (0) 

TO? (</>) = TO^ + TO 2 ((/>) — - Sin 2 COS 2/?TO 2 z( 

Oil ' 3 



TO 



(6) 



and similarly for the fourth generation by substituting with primed quantities, where sin(9yK is 
the weak mixing angle. 

The parameters At, A b , A t i and Ay in the above expressions for the scalar quark masses 
are given as At = At — ix cot (3, and Af, = A b — \x tan (3, and similarly for the fourth generation. 
Note that A t = (Af, = 0) does not necessarily imply that the right-handed and the left-handed 
scalar top (bottom) quarks are degenerate in mass, since there is D-term contributions. Only if 
D-term contributions to the scalar top (bottom) quark masses are neglected, A t = (A b = 0) 
would yield degenerate right-handed and left-handed scalar top (bottom) quarks. 

However, we remark that the parameters A t , A b , A t > and Ay control the mixings between the 
scalar top or scalar bottom masses in each generation. If these parameters are zero, there would 
be no mixing between right-handed and left-handed scalar quarks of each generation. In this 
paper, we assume that there is no mixing, taking A t = Ab = A t > = Ay = in the expressions 
for the scalar quark masses. Therefore, we study the MSSM Higgs potential in the decoupling 
limit without mixing. 

The decoupling limit without mixing is an optimal situation for electroweak phase transi- 
tion to be strongly first order. The case without mixing is more favorable for the first order 
electroweak phase transition to be strong than the case with mixing, as the strength of the 
electroweak phase transition is found to decrease when the mixings between scalar quarks are 
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taken into account [15]. Also, one can notice, for example, in Fig. 2 of Ref. [16], that the 
strength of the electroweak phase transition increases as m\ increases in the MSSM with three 
generations. Thus, we study the MSSM Higgs potential in circumstances that enhance the first 
order electroweak phase transition. 

Now, the renormalized parameter m\ in the Higgs potential can be eliminated by the mini- 
mum condition. By calculating the first derivative of V((f>, 0) with respect to (j), m\ is expressed 
as 



m n = 



cos 2 2(5 + 




1 



(7) 



where a = ti, t 2 , bi, b 2 , and a! = t'i, t' 2 , b'i, b' 2 , and v = 246 GeV. The mass of h at the one-loop 
level in the decoupling limit without mixing at zero temperature is obtained by calculating the 
second derivative of V((f>, 0) with respect to 4> as 



m| cos 2 2(5 + 



3mj 




(8) 



where a = t±, t 2 , bi, b 2 , and a' = t[, t' 2 , b[, b' 2 . 

Now, we study the effect of finite temperature. The one-loop contribution at finite temper- 
ature is given by [17] 



(9) 



V 1 (<f>, T) = J™ dx x 2 log 1 ± exp (-yfx* + m?{(t>)/T* 



where I = W, Z, t, b, ti, t 2 , &i and b 2 , t' , b', t 2 , b'^ and b' 2 . and the negative sign is for bosons 
and the positive sign for fermions. The full one-loop effective potential at finite temperature 
that we are considering can now be expressed as 



V((f>, T) = V (0, 0) + Vi (0, 0) + Vi (</>, T) 



(10) 



We perform the exact integration in V(<p, T) instead of employing the high-temperature approx- 
imation. 
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III. Numerical Analysis 



At the tree level, we have only one free parameter tan (5 in the decoupling limit where tua S> mz- 
At the one-loop level, the number of free parameters increases as m t , nib, mt', my, and thq, 
rriQi are introduced to the one-loop contributions. To be concrete, we set mt = 175 GeV and m& 
= 4.5 GeV, and we take for simplicity the soft SUSY breaking parameters as tuq = mfj = m 2 D , 
and similarly for the fourth generation. Since we assume no mixing in the scalar quark sector, 
we set A t = At, = A t > = Ay = 0. For the masses of the fourth generation quarks, there 
are some experimental constraints. Some years ago, Tevatron data have set my > 119 GeV, 
and recently, from the search for long-lived charged massive particles at Tevatron come more 
stringent experimental lower bounds of my > 180 GeV and m? > 230 GeV [18]. With these 
constraints in mind, we take m t > = 250 GeV and my = 200 GeV. Finally, we take 1 TeV for 
the value of ttiq from the SUSY breaking scale Mstjsy = 1 TeV. Thus, our numerical analysis 
involves two free parameters: tan (5 and mQ> . 

In Fig. 1, we show a typical behavior of V(<j>, T) as a function of <\>, at a critical temperature 
T = T c = 84.225 GeV, where we take tan/3 = 20 and vtlqi = 100 GeV. As one can see in 
the figure, we obtain the critical vacuum expectation value as v c = 149 GeV, and the ratio as 
v c /T c = 1.769. one can notice that the potential in Fig. 1 allows a strongly first-order electroweak 
phase transition. The scalar quark masses are obtained as = 1013 GeV, m ( - 2 = 1014 GeV, 
m- bi = 1000 GeV, = 1001 GeV, mp = 263 GeV, m t ~, = 266 GeV, my = 225 GeV, and 
mr! = 231 GeV. 

For the mass of h, we obtain m^ = 129 GeV with the parameter values of Fig. 1. This 
number is a little larger than the experimental lower bound on the SM Higgs boson mass, 115 
GeV. In the decoupling limit of m\ S> mz, the behavior of h is identical to that of the SM 
Higgs boson. The search for a light Higgs boson in the MSSM by DELPHI collaboration without 
mixing, where only three generations of quarks are taken into account, suggests that mh is about 
115 GeV, for tan/3 > 15 [19]. Comparing this number with our result of 129 GeV in Fig. 1, 
we may well deduce that the difference is due to the contribution by the fourth generation of 
quarks. The contribution by the fourth generation of quarks also enables the electroweak phase 
transition strongly first order. Our value is compatible with the result of Ref. [8] where the 
upper bound on the lightest Higgs boson mass is obtained as 152 GeV for small tan (5 and 
96 < m t i,my < 125 GeV in the MSSM with four generations of quarks. If the mass difference 
between if and b' is small, they may have smaller masses [9]. 

Now, let us study other regions of the parameter space. We vary mni while fixing tan f3 = 20. 
The masses of the third generation of scalar quarks are then the same as Fig. 1: m^ = 1013 GeV, 
mi 2 = 1014 GeV, = 1000 GeV, and = 1001 GeV. The masses of the fourth generation 
of scalar quarks, as well as other relevant quantities, are shown in Table I. The second row 
of Table I for mni = 100 GeV corresponds to the numerical result of Fig. 1. One can see 
that the strength of the first order electroweak phase transition decreases as mQi increases until 
mgt reaches 140 GeV, beyond which v c /T c becomes less than 1. Thus, the electroweak phase 
transition remains strongly first order for mni < 140 GeV. 

We need to explore the boundary of the parameter space beyond which the electroweak phase 
transition is no longer strongly first order. In Table I, one can see that v c /T c ~ 1 for tan/? = 20 
and mQi = 140 GeV. From this point, we examine several points of (mQ>, tan/3) which yield 
v c /T c ~ 1, by adjusting T c . For each value of tan/3, we find the upper bound value of mg/, 
beyond which v c /T c becomes less than 1. The result is shown in Table II. The fourth row of 
Table II is corresponds to the fourth row of Table I. The numbers in Table II indicate that the 
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electroweak phase transition can be strongly first order for 2 < tan f3 < 40 if mn> < 140 GeV. 

In Fig. 2, we plot the numerical results of Table 2 on {mni, tan /3)-plane. The dashed curves 
denote the contours of m^. The solid curve denotes the contour of v c /T c = 1. On the left- 
hand side of the solid curve we have v c /T c > 1, that is, the electroweak phase transition is 
strongly first order in the region to the left of the solid curve. Consequently, we find a region in 
the parameter space of the MSSM with four generations of quarks where a strongly first order 
electroweak phase transition is allowed. The parameter values of the allowed region are within 
the experimental constraints, and yield mh consistent with experimental lower bound. 

Table I: Some values of mqi for which the first order electroweak phase transition is strong 
(v c /T c > 1.0). The critical temperatures T c are obtained for which the finite temperature 
effective potential has two degenerate vacua. The relevant parameter values are the same as 
Fig. 1, that is, tan/3 = 20, m Q = 1 TeV, m t > = 250 GeV, and m h < = 200 GeV. 
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249 
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2.431 


100 


263 
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231 


129 


84.225 


149 


1.769 


130 


276 


279 


239 


245 
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90.180 


113 


1.253 


140 


281 


284 


245 


251 


136 


92.268 


93 


1.007 


150 


286 


289 


251 


256 


138 


94.384 


76 


0.805 



Table II: Some values of (rriQ> , tan (3) beyond which the first order electroweak phase transition 
is strong (v c /T c ~ 1.0). The relevant parameters values are the same as Table I. The masses of 
the scalar quarks of the fourth generation are calculated as about mp = 282 GeV, mp = 285 
GeV, mr, = 246 GeV, my = 251 GeV. 

°1 2 



tan/3 


m Q' 


m h 


T c 


V c 


Vc/Tc 


2 


147 


120 


87.163 


88 


1.009 


5 


142 


133 


91.256 


92 


1.008 


10 


140 


135 


91.965 


92 


1.000 


20 


140 


136 


92.268 


93 


1.007 


30 


140 


137 


92.326 


93 


1.007 


40 


140 


140 


92.346 


93 


1.007 



IV. Conclusions 

Up to now, we study the possibility of a strongly first order electroweak phase transition in the 
MSSM with sequential four generations of quarks. We assume that m^ 3> mz and A\ = 0(7 = 
i, b, t' , b'), that is, we work in the decoupling limit and in the case of no mixing between scalar 
quarks. We choose the relevant parameter values to be: m t = 175 GeV, m^ = 4.5 GeV, m t > = 
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250 GeV, my = 200 GeV, and iuq = 1 TeV. We take mq = mfj = m 2 D , and similarly for the 
fourth generation. These numbers are consistent with experimental constraints. 

We search the parameter space of (mQ/,tan/3)-plane to examine if the electroweak phase 
transition is strongly first order. We find that there are regions in the (mQ> , tan /3)-plane that 
satisfy our criterion of v c /T c > 1. For 2 < tan/3 < 40, the electroweak phase transition is 
strongly first order in the region where rriQi < 140 GeV. The scalar quark masses of the fourth 
generation are controlled mainly by the soft SUSY breaking parameter mq' . In the region where 
the electroweak phase transition is strongly first order, the scalar quark masses of the fourth 
generation are obtained to be larger than the quark masses of the fourth generation. 
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Figure Captions 

Fig. 1: The plot of V(<fi, T) as a function of 4>, for the critical temperature T c = 84.225 GeV. 
The relevant parameter values are set as follows: tan (3 = 20, mg = 1 TeV, mp = 250 GeV, 
my = 200 GeV, and mqi = 100 GeV. One can see that the electroweak phase transition is 
first order. The critical VEV is obtained as v c = 149 GeV. Thus, v c /T c = 1.769, and the first 
order electroweak phase transition is strong. The potential in this figure yields nih = 129 GeV, 
m k = 1013 GeV, m t ~ 2 = 1014 GeV, m~ bi = 1000 GeV, = 1001 GeV, m ? = 263 GeV, 
mp^ = 266 GeV, my = 225 GeV, and m^= 231 GeV. 

Fig. 2: Plots of the lightest scalar Higgs boson mass and the criterion of the strongly first order 
electroweak phase transition in the (mQi, tan /?)-plane. The remaining relevant parameters are 
set as ttiq = 1 TeV, m t > = 250 GeV, and my = 200 GeV. The solid curve is the contour of 
v c /T c = 1, and the dashed curves are the contours of m^ = 115, 120, 125, 130, 135, and 140 
GeV. The masses of the scalar quarks of the third generation are calculated the same as Fig. 1 , 
that is, m h = 1013 GeV, m~ t2 = 1014 GeV, m~ bi = 1000 GeV, and m~ b2 = 1001 GeV. The region 
to the left-hand side of the solid curve and above the dashed curve of mh = 115 GeV is where 
the first order electroweak phase transition is strong (v c /T c > 1) and m^ is consistent with the 
experimental lower bound. 
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Fig. 1: The plot of V((ft,T) as a function of (ft, for the critical temperature T c = 84.225 GeV. 
The relevant parameter values are set as follows: tan/3 = 20, mg = 1 TeV, m t i = 250 GeV, 
my = 200 GeV, and mgi = 100 GeV. One can see that the electroweak phase transition is 
first order. The critical VEV is obtained as v c = 149 GeV. Thus, v c /T c = 1.769, and the first 
order electroweak phase transition is strong. The potential in this figure yields rrih = 129 GeV, 
m fl = 1013 GeV, m t ~ 2 = 1014 GeV, m~ bi = 1000 GeV, = 1001 GeV, m t ~, = 263 GeV, 
m t ~, = 266 GeV, my = 225 GeV, and m^= 231 GeV. 
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Fig. 2: Plots of the lightest scalar Higgs boson mass and the criterion of the strongly first order 
electroweak phase transition in the (rriQi, tan /3)-plane. The remaining relevant parameters are 
set as wiq = 1 TeV, m t i = 250 GeV, and my = 200 GeV. The solid curve is the contour of 
v c /T c = 1, and the dashed curves are the contours of rrih = 115, 120, 125, 130, 135, and 140 
GeV. The masses of the scalar quarks of the third generation are calculated the same as Fig. 1, 
that is, m h = 1013 GeV, m fa = 1014 GeV, m~ bi = 1000 GeV, and = 1001 GeV. The region 
to the left-hand side of the solid curve and above the dashed curve of m-h = 115 GeV is where 
the first order electroweak phase transition is strong (v c /T c > 1) and rrih is consistent with the 
experimental lower bound. 
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